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PREFACE 


This  report  describes  a method  for  estimating  wind-wave  growth  and 
decay  over  flooded  areas  where  there  is  a major  friction  effect  because 
of  dense  vegetation.  The  report  was  initiated  in  response  to  a request 
from  the  U.S.  Army  Engineer  Division,  Lower  Mississippi  Valley,  New 
Orleans  District  at  the  Division’s  14  September  1976  Research  and  Devel- 
opment Workshop,  indicating  a need  for  technical  guidelines  for  predicting 
wind-wave  generation  over  flooded  coastal  areas.  The  work  was  carried 
out  under  the  coastal  construction  program  of  the  U.S.  Army  Coastal 
Engineering  Research  Center  (CERC) . 

These  technical  guidelines  are  an  extension  of  the  procedures  given 
in  the  Shore  Protection  Manual  (SPM)  (U.S.  Army,  Corps  of  Engineers, 
Coastal  Engineering  Research  Center,  1975).  The  design  curves  in  the 
SPM  are  limited  to  waves  passing  over  a sandy  bottom.  A condensed  des- 
cription of  the  method  is  presented  in  CERC  Coastal  Engineering  Technical 
Aid  No.  77-6  (Camfield,  1977). 

This  report  was  prepared  by  Dr.  Frederick  E,  Camfield,  Hydraulic 
Engineer,  under  the  general  supervision  of  R.A.  Jachowski,  Chief,  Coastal 
Design  Criteria  Branch. 

Comments  on  this  publication  are  invited. 


Approved  for  publication  in  accordance  with  Public  Law  166,  70th 
Congress,  approved  31  July  1945,  as  supplemented  by  Public  Law  172, 
88th  Congress,  approved  7 November  1963. 
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CONVERSION  FACTORS,  U.  S,  CUSTOMARY  TO  METRIC  (SI) 
UNITS  OF  MEASUREMENT 


U.S.  customary  units  of  measurement  used  in  this  report  can  be  converted  to  metric  (SI) 
units  as  follows: 


To  obtain 


inches 

25.4 

millimeters 

2.54 

centimeters 

square  inches 

6.452 

square  centimeters 

cubic  inches 

16.39 

cubic  centimeters 

feet 

30.48 

centimeters 

0.3048 

meters 

square  feet 

0.0929 

square  meters 

cubic  feet 

0.0283 

cubic  meters 

yards 

0.9144 

meters 

square  yards 

0.836 

square  meters 

cubic  yards 

0.7646 

cubic  meters 

miles 

1.6093 

kilometers 

square  miles 

259.0 

hectares 

knots 

1.8532 

kilometers  per  hour 

acres 

0.4047 

hectares 

foot-pounds 

1.3558 

newton  meters 

millibars 

1.0197  X 10'^ 

kilograms  per  square  centimeter 

ounces 

28.35 

grams 

pounds 

453.6 

grams 

0.4536 

kilograms 

ton,  long 

1.0160 

metric  tons 

ton,  short 

0.9072 

metric  tons 

degrees  (angle) 

0.1745 

radians 

Fahrenheit  degrees 

5/9 

Celsius  degrees  or  Kelvins' 

To  obtain  (.rigiu.s  ((.)  temperature 
To  obtain  Kelvin  (K)  readines,  use 


readinfip  from  rahrenheit  (r)  readiiipt,  use  formula:  - (o/O)  (F  - ,‘$2) 
formula:  K = (5/9)  (K  - 32)  + 273.1 5. 


SYMBOLS  AND  DEFINITIONS 


water  depth 

water  depth  at  seaward  or  beginning  edge  of  segment 
fetch  length 

adjusted  fetch  length  for  distance  across  a segment  in  the 
direction  of  wave  motion 

equivalent  fetch  length  for  the  initial  wave  at  seaward  or 
beginning  edge  of  segment 

bottom-friction  factor  (Darcy-Weisbach  friction  factor) 
bottom-friction  factor  at  seaward  or  beginning  edge  of  segment 
fractional  growth  factor  of  equivalent  initial  wave 
gravitational  acceleration 
wave  height 

decayed  wave  height  at  end  of  fetch 
equivalent  wave  height  at  end  of  fetch 
wave  height  at  end  of  fetch 

wave  height  at  seaward  or  beginning  edge  of  segment 
equivalent  initial  wave  height 
maximum  stable  wave  height 

maximum  significant  wave  height  which  would  be  generated  for 
given  windspeed  and  water  depth 

decay  factor 

decay  factor  when  the  bottom-friction  factor,  f^  = 0.01 

decay  factor  when  the  bottom-friction  factor,  ff,  has  a value 
different  than  0.01 

shoaling  coefficient 

wavelength 

deepwater  wavelength  for  a wave  with  period,  T 


SYMBOLS  AND  DEFINITIONS--Continued 


n Manning's  roughness  coefficient 

p incremental  multiple  of  wave  height 

fractional  reduction  of  initial  wave  at  seaward  edge  of  segment, 
as  compared  to  the  maximum  stable  wave  height 

T wave  period 

U windspeed 

X distance  in  the  direction  of  wave  motion 

a factor  for  reducing  fetch  length  to  the  adjusted  length 

Op  factor  for  increasing  fetch  length  to  the  adj listed  length  = l/o 

A incremental  change 

Ax  actual  distance  across  a segment  in  the  direction  of  wave  travel 


WIND-WAVE  PROPAGATION  OVER  FLOODED,  VEGETATED  LAND 

Frederick  E.  Ccmt field 
I . INTRODUCTION 

An  important  factor  in  the  planning  and  design  of  works  to  protect 
upland  property  during  periods  of  storm  surge  is  the  prediction  of  the 
wave  height  and  period  that  will  prevail  at  and  seaward  of  the  protec- 
tive works  (i.e.,  levee,  dike,  seawall)  for  a selected  design  storm. 
Although  improvements  are  needed,  guidelines  are  available  for  predic- 
tion of  the  water  levels  in  upland  areas  resulting  from  storm  surge; 
however,  no  guidelines  are  presently  available  for  computing  the  wave 
attenuation  for  conditions  when  storm- generated  waves  travel  a distance 
across  shallow  flooded  areas  where  the  bottom  characteristics  include 
vegetation  which  causes  a moderate  to  high  frictional  stress.  Therefore, 
it  is  necessary  to  estimate  the  heights  and  periods  of  these  storm- 
generated waves. 

Theoretical  methods  and  field  measurements  are  currently  unavailable 
to  provide  an  estimation  procedure.  This  report  presents  a prelimincay 
(approximate)  method  for  estimating  the  growth  or  decay  of  waves  travel- 
ing through  shallow  water  over  areas  with  a high  frictional  resistance 
from  vegetation.  The  method  is  based  on  previously  developed  equations 
(U.S.  Army,  Corps  of  Engineers,  Coastal  Engineering  Research  Center,  1975) 
for  wave  growth  over  areas  with  low  bottom  friction,  and  an  equation  for 
the  decay  of  gravity  waves  over  areas  with  a constant  water  depth  and 
high  bottom  friction.  This  method  uses  existing  shallow-water  wave  fore- 
casting curves  by  adjusting  fetch  lengths  to  account  for  higher  bottom 
friction.  Simplifying  assumptions  are  used.  The  water  depth  is  assumed 
to  have  only  gradual  variations,  and  the  frictional  resistance  is  treated 
as  bottom  friction.  The  procedure  discussed  in  this  report  has  not  been 
verified  in  the  field  and  may  not  be  applicable  to  other  problems  relat- 
ing frictional  resistance  to  wave  development. 

Only  limited  data  are  available  on  the  effects  of  high  values  of 
bottom  friction  on  wind  waves.  Friction  factors  are  estimated  by  compar- 
ing vegetation  to  similar  conditions  in  river  channels  and  on  flood  plains. 
The  effect  of  the  vegetation  on  wind  stress  and  the  possible  effects  of 
motion  of  the  vegetation  are  not  considered.  Dense  vegetation  effects 
near  the  water  surface  which  will  dampen  short-period  waves  much  faster 
than  long-period  waves  are  also  not  considered.  The  results  obtained 
are  considered  to  be  conservative;  i.e.,  the  predicted  wave  heights  are 
expected  to  be  slightly  higher  than  the  wave  heights  which  actually  occur. 

II.  WAVE  FORECASTING  CURVES 

The  Shore  Protection  Manual  (SPM)(U.S.  Army,  Coastal  Engineering 
Research  Center,  1975)  discusses  a procedure  for  predicting  the  growth 
of  waves  over  shallow  water,  which  is  a modification  of  a method 


9 


developed  by  Bretschneider  (1952,  1958,  1970)  and  incorporates  the 
results  of  Ijima  and  Tang  (1966).  The  SPM  method  assumes  that  bottom 
friction  ff  - 0.01.  The  significant  wave  height  and  wave  period  are 
obtained  from  empirical  equations  given  as 


and 


where 

H = significant  wave  height 
T = significant  wave  period 
d = water  depth 
U = windspeed 
F = fetch  length 
g = gravitational  acceleration. 

Since  the  terms  are  evaluated  in  dimensionless  form,  any  consistent  sys- 
tem of  measurements  may  be  used.  Equations  (1)  and  (2)  are  plotted  as 
curves  in  Figures  1 and  2,  respectively,  for  constant  values  of  water 
depth.  Design  curves  for  predicting  the  growth  of  waves  in  shallow  water 
for  specified  constant  water  depths  are  shown  in  Figures  3 to  12  (after 
SPM,  U.S.  Army,  Corps  of  Engineers,  Coastal  Engineering  Research  Center, 
1975) . 

III.  CALCULATION  OF  ADJUSTED  FETCH 

The  long  dashline  in  Figure  1 indicates  the  maximum  significant  wave 
heig^it,  which  will  occur  for  any  given  water  depth  and  windspeed 
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= 10  Osec 


when  the  fetch  is  sufficiently  long.  From  Figure  1,  the  relationship 
between  the  wave  height,  and  the  water  depth,  d,  is  given  by 

iJkil=  0.146  ( 


Where  the  bottom- friction  factor  fy>  = 0.01,  and  for  a given  water  depth 
and  windspeed,  when  the  initial  wave  hei^t,  H^,  at  the  beginning  of 
the  fetch  is  less  than  Hg^,  the  wave  will  continue  to  grow  in  height  as 
it  travels  across  the  fetch,  with  the  height  approaching  the  maximum  sig- 
nificant wave  height,  Hg;^.  When  is  greater  than  Hg;„,  it  is  as- 

sumed that  the  wave  hei^t  will  decay  as  the  wave  travels  across  the 
fetch,  with  the  wave  height  again  approaching  the  maximum  significant 
wave  height,  Hg^. 

Where  the  bottom- friction  factor  ff  > 0.01,  i.e.,  vegetation  in- 

creases the  frictional  resistance,  and  the  initial  wave  heigjit 
it  is  assumed  that  the  wave  will  grow  to  a final  height  less  than  or 
equal  to  Hg^,  but  that  the  growth  will  occur  at  a slower  rate  than 
where  ~ 0*01  (Fig-  13).  Althou^  the  hi^er  frictional  resistance 
may  limit  the  wave  growth  in  this  case  to  a height  less  than  Hg^,  no 
data  are  available  on  this  effect.  Therefore,  as  a conservative  estimate, 
it  is  assumed  that  the  maximtun  significant  wave  height  will  have  the  same 


value, 


as  where  f^  = 0.01. 


Where  the  bottom- friction  factor  f f > 0.01,  and  the  initial  wave 
height  > Hg^,  it  is  assumed  that  the  wave  will  decay  to  a final 
heigjit  less  than  or  equal  to  Hg^,  but  that  the  decay  will  occur  at  a 
greater  rate  than  where  0.01  (Fig.  14).  Althou^  the  hi^er  fric- 

tional resistance  may  cause  the  wave  to  decay  to  a height  less  than  ttgrrit 
no  data  are  presently  available.  Therefore,  the  maximum  significant  wave 
height  is  assumed  equal  to  the  value  of  Hg^  where  ff  = 0.01  as  before. 


Where  the  bottom- friction  factor 
justed  to  an  equivalent  fetch,  F^, 


fj?  > 0.01,  the  fetch  length  is  ad- 
which  has  a bottom  friction  ff  = 


0.01,  and  which  will  give  the  same  growth  or  decay  as  the  actual  fetch 
with  the  actual  bottom  friction.  The  relationship  between  the  adjusted 


fetch. 


and  the  actual  fetch  distance.  Ax,  is  shown  in  Figures  13 


and  14.  After  the  adjusted  fetch,  F^,  has  been  determined.  Figures  1 to 
12  are  used  to  predict  wave  growth  or  decay  (see  Secs.  IV  and  V). 

1.  Determination  of  Friction  Factor. 

Limited  data  are  available  to  define  friction  factors  for  water  wave 
motion  through  dense  stands  of  marsh  grass,  brush,  or  trees.  Saville 
(1952)  presented  marsh  correction  factors  for  correcting  values  of  wave 
setup  over  flooded  marsh  areas.  However,  this  was  not  extended  to  pro- 
vide corrections  for  wave  heights.  Whitaker,  et  al.  (1975)  developed 
numerical  simulations  of  water  level  changes  which  considered  the  effects 


of  vegetation  in  terras  of  drag  coefficients,  but  effects  on  wave  heights 
were  not  considered.  Wayne  (1975)  investigated  the  decay  of  very  low 
aji^Jlitude  waves,  from  0.126  to  0.744  foot  high  (0.038  to  0.227  meter), 
traveling  short  distances  over  grass.  However,  he  did  not  establish  the 
dependence  of  wave  period,  initial  wave  height,  and  water  depth,  or  con- 
sider the  corabination  of  wind  stress  with  bottom  friction. 

To  obtain  values  of  bottom- friction  factors,  estimates  of  Manning's 
roughness  coefficient,  n,  have  been  made  using  values  shown  in  Chow 
(1959)  for  flow  over  flood  plains.  These  estimated  roughness  coefficients 
are  conservative;  i.e.,  they  are  expected  to  give  predicted  wave  heights 
somewhat  greater  than  the  wave  heights  which  actually  occur.  The  re- 
duction in  wind  stress  due  to  vegetation  extending  above  the  water  sur- 
face is  not  considered.  The  roughness  coefficient,  n,  can  be  related 
to  the  friction  factor,  fj>,  for  flow  over  ground  by  combining  the 
Darcy-Weisbach  and  Manning  equations  for  energy  losses  due  to  friction. 
However,  the  Manning  equation  contains  a dimensional  coefficient  which 
varies  depending  on  the  dimensional  units  in  the  equation.  Using  a foot- 
pound-second system  of  units. 


f^  = 


3.60 


vF- 


(4) 


where  d is  in  feet  and  g is  in  feet  per  second  squared.  (In  a meter- 
kilogram-second  system  of  units,  the  coefficient  3.60  in  equation  (4)  be- 
comes 8.0.) 

The  roughness  coefficient,  n,  will  vary  as  a function  of  depth  due 
to  the  height  of  the  vegetation  in  relation  to  the  depth  of  the  water. 

In  addition,  the  friction  factor,  ff,  has  an  inverse  relationship  to 
the  water  depth,  d,  as  shown  by  equation  (4).  Values  of  f«  used  for 
various  kinds  of  ground  cover  are  shown  as  curves  A to  D in"^  Figure  15. 
Curve  A is  for  a sandy  bottom  (Bretschneider , 1952,  1958,  1970),  where 
ff  is  assumed  constant.  This  curve  was  used  for  developing  the  curves 
in  Figures  1 to  12.  Curve  B is  for  coastal  areas  with  thick  stands  of 
grass.  The  grass  will  create  a high  resistance  at  low  water  levels  where 
the  depth  of  water  and  the  height  of  the  grass  are  nearly  equal,  but  the 
resistance  will  rapidly  diminish  as  the  water  depth  becomes  greater  than 
the  grass  height;  i.e.,  the  water  particle  motion  under  the  wave  is  in- 
fluenced less  by  the  vegetation.  Curve  C is  for  higher  levels  of  veg- 
etation, such  as  brush  or  low  bushy  trees  which  extend  above  the  hei^t 
of  the  grass.  The  friction  effects  will  be  higher  for  curve  C than  for 
curve  B because  the  hi^er  vegetation  will  have  more  influence  on  the 
water  particle  motion  for  any  given  water  depth.  Curve  D is  for  tall 
trees  which  will  always  be  higher  than  the  depth  of  water.  The  frictional 
resistance  will  be  high  at  all  water  levels,  and  somewhat  higher  at  low 
water  levels  because  of  the  added  resistance  from  the  ground.  Curve  D 
is  for  a relatively  close  spacing  of  trees;  e.g.,  a second-growth  pine 
forest  which  will  give  the  hipest  frictional  resistance.  A more  scat- 
tered spacing  of  trees  would  give  a lower  value  of  ff  for  a given 


Friction  Foctor  (fj) 


water  depth.  Chow  (1959)  provides  values  of  n for  a wider  range  of 
conditions. 

2.  Selection  of  Fetch  Segment. 

The  growth  or  decay  of  a wave  at  any  point  is  dependent  on  the  water 
depth,  the  wave  height  and  period  at  that  point,  the  bottom  friction,  and 
the  windspeed.  For  the  method  used  here,  the  windspeed  is  assumed  to  be 
constant  across  a fetch.  To  accurately  predict  the  growth  or  decay  of 
waves  traversing  a particular  fetch,  it  is  necessary  to  divide  the  fetch 
into  segments  according  to  water  depth  and  bottom  friction.  The  bottom 
friction  may  vary  substantially  as  a function  of  water  depth  for  a par- 
ticular type  of  vegetation  (Fig.  15). 

Figures  1 to  12  were  derived  for  a constant  depth;  any  variation  in 
depth  is  assumed  to  be  very  gradual,  and  these  figures  are  applied  to  an 
average  depth  across  a fetch  segment.  Therefore,  the  total  variation  in 
depth  across  a fetch  must  be  considered.  The  type  of  vegetation  may  also 
vary,  with  sections  of  marsh  grass,  brush,  trees,  or  shallow  lagoons. 

Wave  heists  will  normally  vary  across  a fetch  and,  as  the  decay  factors 
will  depend  on  the  wave  heig^it,  new  decay  factors  must  be  calculated  if 
the  wave  height  varies  excessively.  A fetch  segment  is  also  considered 
to  be  much  longer  than  a single  wavelength. 

Dividing  the  fetch  into  segments,  the  segment  distance  (length  in  the 
direction  of  wave  travel)  is  determined  so  that,  first. 

Ad  < 0.25  di  , (5) 

where  Ad  is  the  change  in  depth  over  the  segment  distance  and  d^  is 
the  depth  at  the  seaward  or  beginning  edge  of  the  segment;  second, 

Af^  < 0.25  f^^  , (6) 

where  Afy?  is  the  change  in  the  bottom- friction  factor  over  the  seg- 
ment distance  and  fy^  is  the  friction  factor  at  the  seaward  or  begin- 
ning edge  of  the  segment;  and  third,  after  the  change  in  wave  height 
across  the  segment,  AH,  has  been  determined, 

AH  < 0.5  , (7) 

where  H^  is  the  initial  wave  hei^t  at  the  seaward  edge  of  the  segment. 
Wave  decay  is  also  assumed  to  be  small  (less  than  0.1  H^)  over  a single 
wavelength.  The  segment  distance  may  require  adjustment;  i.e.,  a shorter 
segment  may  be  necessary  if  the  confuted  value  of  AH  is  unacceptable. 
The  numerical  coefficients  used  on  the  rigjit-hand  side  of  equations  (5)  , 
(6),  and  (7)  are  arbitrarily  chosen.  Smaller  coefficients  would  be  ex- 
pected to  increase  the  calculations  required  for  a solution;  larger  co- 
efficients would  be  expected  to  produce  a greater  error  in  the  estimated 
wave  height  obtained. 


Because  of  the  storm  conditions  producing  the  surge,  it  is  assumed 
that  waves  will  always  exist  at  the  seaward  edge  of  the  initial  segment 
(i.e,,  >>  0),  and  values  of  wind  velocity,  U,  greater  than,  e.g., 

50  miles  per  hour  (73  feet  or  22.4  meters  per  second)  will  be  expected. 

In  case  where  = 0,  the  methods  presented  here  do  not  apply.  Initial 
wave  growth  would  consist  of  very  low  amplitude  waves  with  very  short 
periods.  These  initial  waves  will  only  be  influenced  by  vegetation  if 
the  vegetation  extends  above  the  wave  surface.  That  case  has  not  been 
investigated,  and  cannot  be  treated  by  the  method  discussed  in  this 
report . 

3.  Adjustment  of  Fetch  Length. 

The  adjusted  fetch  length  of  the  segment  is  determined  by  comparing 
the  decay  factor,  Kf,  for  the  actual  bottom  friction  of  the  segment 
(from  Fig.  15)  with  the  decay  factor  obtained  for  the  friction  factor 
ff  = 0.01  which  was  used  for  the  wave  predictions  in  Figures  1 to  12. 
Using  the  method  of  Bretschneider  (1954)  (Fig.  16),  the  decay  factor, 

Kf,  is  obtained  by  using  the  two  factors  f^  Ax/d^  and  2iTd/(gT^), 
where  Ax  is  the  actual  segment  distance  in  feet  in  the  direction  of 
wave  motion. 

The  decay  factor,  K^,  was  defined  by  Bretschneider  as 


ff  Hi  Ax 

/ 27rd 

16TT  K5.2  , . 

d2 

7 

(sinh  ’ 

(8) 


where  Kg  is  the  shoaling  coefficient  defined  as  in  Table  C-1  of 

the  SPM.  By  noting  that 


(9) 


where  is  the  deepwater  wavelength,  and  Kg  and  sinh  (2iTd/L) 

are  functions  of  d/L^  and  therefore  of  2Trd/(gT2),  it  can  be  shown 
that  when  2iTd/(gT^)  is  held  constant  that  the  right-hand  side  of  equa- 
tion (8)  is  then  a function  of  f^  Ax/d^  only. 


Bretschneider  (1954)  used  the  work  of  Putnam  and  Johnson  (1949)  to 
define  a friction  factor  for  an  impermeable  sandy  bottom  as  ~ 0.01. 
The  curves  in  Figures  1 to  12  were  developed  for  that  friction  factor. 
The  wave  decay  is  proportional  to  (1  - Kj?) . To  compare  wave  decay  over 
a fetch  segment  with  high  bottom  friction  to  wave  decay  over  a fetch 
segment  where  the  bottom  friction  f^  = 0.01  (\ii,  d,  and  U being  the 

same),  a proportionality  factor,  a,  is  defined  by 
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a 


(10) 


(1  - ^f.oi) 
(1  - K/a) 


where  K^.  qj  is  the  decay  factor  defined  for  f^  = 0.01  and  is  the 

decay  factor  defined  for  the  actual  bottom  friction  with  water  depth,  d, 
initial  wave  height,  H^,  windspeed,  U,  and  segment  length.  Ax,  re- 
maining the  same. 

For  wave  growth,  a segment  with  bottom  friction  f^  > 0.01  would 
require  a longer  segment  length  for  waves  to  grow  to  a given  height  than 
the  length  that  would  be  required  where  f^  = 0.01.  Therefore,  the  seg- 
ment with  greater  bottom  friction  is  equivalent  to  a shorter  segment  with 
fy  = 0.01.  As  an  approximation,  the  shorter  adjusted  segment  length,  Fq, 
is  defined  as 


Pa  = 


a Ax 


(11) 


where  Ax  is  the  actual  segment  length  (see  Fig.  13) . 


For  wave  decay,  a fetch  segment  with  bottom  friction  fy  > 0.01 
would  require  a shorter  length  for  waves  to  decay  to  a given  height  than 
would  be  required  where  f^  = 0.01.  Therefore,  the  segment  with  higher 
bottom  friction  is  equivalent  to  a longer  segment  with  f^  = 0.01  (see 
Fig.  14).  The  longer  adjusted  segment  length,  F^,  in  this  case  is 
defined  as 


F^  = Uj,  Ax  , 


where 


1 

a - (1  - Ky.o]) 


(12) 


(13) 


IV.  WAVE  GROWTH  IN  SHA' LOW  WATER 

From  Figure  1 or  equation  (3),  for  any  given  water  depth,  windspeed, 
and  fetch  length,  a maximum  significant  wave  height.  which  would 

be  generated  can  be  defined.  If  the  initial  wave  height,  H^,  at  the 
seaward  or  beginning  edge  of  the  fetch  segment  is  less  than  Hg^,  it  is 
assumed  that  the  wave  will  grow  to  a higher  height  as  discussed  previously. 

To  determine  the  wave  growth,  it  is  necessary  to  first  determine  an 
equivalent  fetch  length,  Fg,  for  the  initial  wave.  This  is  obtained 
directly  from  Figures  1 to  12  using  the  given  windspeed  and  water  depth. 
Secondly,  the  adjusted  fetch,  F^j,  is  determined  using  equations  (10) 
and  (11)  and  Figure  16.  The  total  fetch  is  then  given  as 


i 


J 


I 


Re-entering  Figures  1 to  12  with  the  fetch  length,  F,  the  windspeed, 
U,  and  water  depth,  d,  the  final  wave  height  at  the  end  of  the  fetch 
segment,  Hj?,  is  determined.  This  is  shown  schematically  in  Figure  17 
and  in  the  following  example  problem. 
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EXAMPLE  PROBLEM  i************ 


GIVEN ; A wave  passes  into  shallow  water  over  a flooded  coastal  area. 

The  water  depth,  d^,  at  the  seaward  edge  of  the  area  is  23  feet  (7 
meters),  and  at  the  landward  edge  of  the  area  the  depth  is  13  feet  (4 
meters).  The  distance  across  the  area  in  the  direction  of  wave  motion 
is  10,000  feet  (3,050  meters).  The  wave  height,  H^,  at  the  seaward 
edge  of  the  area  is  limited  by  large  sandbars  seaward  of  the  area  being 
considered  and  is  3 feet  (0.91  meter);  the  wave  period  is  3.2  seconds. 
The  windspeed  is  70  miles  per  hour  (102.7  feet  or  31.3  meters  per  second) 
The  flooded  area  is  covered  with  thick  stands  of  tall  grass. 

FIND:  The  height  and  period  of  the  significant  wave  at  the  landward  edge 
o?  the  segment. 


SOLUTION: 


0.25  <Li  = 0.25  (23)  = 5.75  feet 


Ad  = 23  - 13  = 10  feet  >0.25  d^  . 

Since  this  does  not  meet  the  condition  of  equation  (5) , the  area  should 
be  divided  into  two  fetch  segments.  Assuming  a uniform  variation  in 
depth,  take  the  first  segment  as  a distance  Ax  = 5,000  feet  with  a 
depth  variation  from  23  to  18  feet.  Then 

Ad  = 23  - 18  = 5 feet  < 0.25  d^  . 

At  the  23-foot  depth  (from  Fig.  15,  curve  B) , 

= 0.080 

and  at  the  18- foot  depth  (curve  B) , 
ff  = 0.095 

Afy  = 0.095  - 0.080  = 0.015 
0.25  ffi'=  0.25  (0.080)  = 0.020 
= 0.25  . 

Equations  (5)  and  (6)  are  satisfied,  so  the  fetch  segment  chosen  is  used. 
For  a uniformly  varying  depth,  the  average  depth  can  be  taken  as  the 
average  of  the  depths  at  the  beginning  and  the  end  of  the  segment;  i.e.. 


; 

i 


j 

I 

( 


d 


23  f 18 
2 


20.5  feet  . 


For  a uniform  type  of  vegetation,  the  friction  factor  will  vary  as  a 
function  of  water  depth  (Fig.  IS) . As  an  approximation,  the  average 
friction  factor  can  be  taken  as  the  average  of  the  friction  factors  at 
the  beginning  and  the  end  of  the  segment,  i.e., 


0.080  * 0.095  ^ 0.088 


For  d » 20.5  feet,  = 3 feet,  and  U * 70  miles  per  hour  from 
Figures  1 or  6 


32.2  X 20.5 
(102.7)2 


0.0626 


8Hi 

U2 


32.2  X 3 

= 0.00916 

(102.7)2 


and  from  Figure  1 


For  ff  = 0.01, 


32.2 


000  feet 


„ Q.Ol  X 3 X 5,000  , 
d2  20.52 


for  fy»  » 0.088, 

H 0.088  X 3 X 5,000  , - . 

d2  “ 20.52 

For  the  period,  t = 3.2  seconds,  and  d > 20.5  feet, 

0.391  . 

gfr2  32.2  (3.2)2 

For  2itd/(gT2)  ■ 0.391  (from  Fig.  16) 
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= 0.96S  for  fj?  = 0.088  and  fy  H£  Ax/d^  = 3.14  . 


From  equation  (10), 


1 - K, 


'.01  1 - 0.996  0.004  „ 

fa  ' 1 - 0.965  0.035  ' ' 


from  equation  (11), 


Fa  = o Ax  = 0. 114  (5,000)  = 570  feet;  and 
from  equation  (14), 

**  = + Pa  = + 570  = 4,570  feet  . 

r 

J 

For  d = 20.5  feet,  U = 70  miles  per  hour,  and  F = 4,570  feet  (from 
Figs.  1 and  2,  or  6) 

Hy  = 3.17  feet  and  T = 3.31  seconds 

AH  = 3.17  - 3 = 0.17  feet  < 0.50  . 

This  satisfies  the  requirements  of  equation  (7),  and  the  solution  may 
proceed  to  the  next  segment  which  is  the  remaining  5,000  feet  of  the 
area,  with  the  water  depth  varying  from  18  to  13  feet,  so 

0.25  d^  = 0.25  (18)  = 4.50  feet  . 

Since  Ad  = 18  - 13  = 5 feet  > 0.25  d^,  which  does  not  satisfy  equation 
(5) , a shorter  segment  is  required.  For  a 3,000- foot  segment,  assuming 
a uniform  depth  variation,  the  depth  will  vary  from  18  to  15  feet.  For 
the  15-foot  depth  (using  curve  B in  Fig.  15) 

fy  = 0.120 

fy^  = 0.095  at  the  18- foot  depth  as  previously  shown. 

Af^  = 0.120  - 0.095  = 0.025  = 0.25  f^v  . 


This  satisfies  equation  (6)  and  the  solution  may  proceed.  The  average 
depth,  d = 16.5  feet,  and  the  average  friction  factor,  fy  = 0.108. 

For  d = 16.5  feet  and  11  = 3.17  feet  (from  Fig.  1) 

F^  = 5,400  feet  ; 


for  d = 16.5  feet,  = 3.17  feet,  f*  = 0.108,  Ax  = 3,000  feet,  and 
T = 3.31  seconds  (from  Fig.  16) 

= 0.294 
8T2 

V-01  = ^ 

= 0.88  for  fy?  = 0.108  and  f^  H/  Ax/d^  = 3.77  .. 

Using  equation  (10),  a = 0.1  and 

F(X  - a Ax  = 0.1  (3,000)  = 300  feet 
F = Fg  + Fa  = 5,400  + 300  = 5,700  feet  . 

For  d = 16.5  feet  (from  Figs.  1 and  2) 

= 3.27  feet  and  T = 3.41  seconds  . 

The  remaining  2,000  feet  of  the  fetch  can  then  be  treated  as  a third 
segment.  The  average  depth,  d = 14  feet,  and  the  average  friction 
factor,  f^  = 0. 13. 

For  d = 14  feet  and  Hi  = 3.27  feet  (from  Fig.  1), 

Fg  = 7,200  feet;  ; 

for  d = 14  feet,  = 3.27  feet,  fy  = 0.13  (from  Fig.  16) 

Ax  = 2,000  feet,  T = 3.41  seconds,  and  2Trd/(gt^)  = 0.235  , 
V-Ol  “ 

Kfa  = 0.80  for  f/  = 0.13  and  ff  Hi  Ax/d^  = 4.34  . 

Using  equation  (10),  a = 0.1  and 

Fg  = a Ax  = 0.1  (2,000)  = 200  feet 

^ = '^e  * ^'200  + 200  = 7,400  feet  . 

For  d = 14  feet,  U = 70  miles  per  hour,  and  F = 7,400  feet  (from  Figs. 
1 and  2) 

Hjf  = 3.34  feet  and  T ® 3.51  seconds  . 

Note. --For  a sandy  bottom  ff  = 0.01,  the  wave  would  have  increased 
to  a height  of  approximately  4.26  feet,  a 42-percent  increase  from 


the  initial  wave  hei^t  of  3 feet.  For  the  thick  stands  of  tall  grass, 
the  predicted  increase  in  wave  heigjit  is  only  11  percent  using  the 
approximate  method  of  solution  discussed  in  this  report. 

**********ilr**il'**il>*ilr****ilr**it>*ilr******.  ** 

V.  WAVE  DECAY  IN  SHALLOW  WATER 

Values  of  wind- generated  significant  wave  heights  and  wave  periods 
as  a function  of  windspeed  and  water  depth  are  shown  in  Figures  1 to  12. 

If  the  initial  significant  wave  height  at  the  seaward  or  beginning  edge 
of  a segment  of  fetch  exceeds  the  maximum  significant  wave  heigjit  for 
the  given  water  depth  of  the  segment  of  fetch  and  the  given  windspeed, 
the  effects  of  the  bottom  friction  may  be  assumed  to  exceed  the  effects 
of  the  wind  stress.  Therefore,  the  wave  is  assumed  to  decay,  lose  height, 
and  over  a long  distance  to  approach  a height  equal  to  the  maximum  sig- 
nificant wave  height. 

The  higgler  waves  at  the  beginning  of  the  fetch  will  actually  repre- 
sent a spectnun  of  waves.  Waves  of  various  periods  could  break  at  dif- 
ferent points  approaching  the  shallow-water  fetch  area  from  deep  water. 

As  shown  in  Figure  16,  the  method  used  here  predicts  that  the  longer 
period  wind  waves  would  decay  faster  than  shorter  period  waves  passing 
throu^  the  shallow  water  (see  Fig.  18).  As  discussed  previously,  the 
effects  of  vegetation  near  the  surface  may  have  a strong  influence  on 
very  short-period  waves.  However,  these  waves  are  considered  to  be  less 
in^ortant  for  design  purposes,  and  this  influence  is  ignored. 


Assumed  Exponential  Decoy  of  Significant  Wove 
_ ^ Decoy  of  Short-Period  Waves 


uecuy  01  onori 
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Fetch  Length, X-  - 

Figure  18,  Exponential  wave  growth  and  decay. 
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As  a first  approximation,  the  wave  decay  may  be  assumed  to  occur 
exponentially  at  a rate  similar  to  the  rate  of  wave  growth;  i.e.,  the 
wave  height  will  decay  from  a maximum  stable  wave  height  given  by  the 
shallow-water  condition 

= 0.78  d (15) 

to  a height  equal  to  the  maximum  significant  height  in  a distance  approx- 
imately equal  to  the  shallow-water  fetch  length  required  for  a wave  to 
grow  from  a zero  height  to  the  maximum  generated  significant  wave  height 
(the  fetch  length  defined  by  the  long  dashlines  in  Fig.  1). 

1.  Wave  Hei^t. 


The  means  of  determining  the  decayed  wave  height  is  shown  schematic- 
ally in  Figure  19.  Steps  to  predict  the  decay  of  a wave  are: 


G i { Hfn”  Hsm  ) 


I u (hi 


0 0 


'3) 

Hj  = Gi Hsm 


Step  in  Colculafions 
Procedure 


Hi*=  Rj  Hjm 


0 Fe 


Figure  19.  Schematic  of  wave  decay  calculations. 


(a)  Determine  the  maximum  significant  wave  height  that  would  be 
generated  for  the  given  windspeed  and  water  depth,  assuming  an  unlimited 
fetch,  and  using  Figures  1 to  12  or  equation  (3). 

(b)  Determine  the  fractional  reduction,  R^,  represented  by  the  initial 
wave  at  the  seaward  edge  of  the  segment  of  fetch  under  consideration 
given  by 


fj*  initm 
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Hm  - ^ 
Hm  - H 


where  is  the  maximum  stable  wave  height  given  by  equation  (15), 

the  incident  wave  height  at  the  seaward  edge  of  the  fetch  segment,  and 
the  maximum  significant  wave  height. 


(c)  Determine  the  equivalent  initial  wave  height,  for  wave 

growth  by 

(d)  Determine  the  equivalent  fetch  length,  Fg,  for  the  wave  height, 

Htg. 


(e)  Determine  an  adjusted  fetch  length,  F^,  for  the  segment  length. 
Ax,  using  equations  (12)  and  (13). 

(f)  Determine  the  total  fetch,  F,  from  equation  (14). 

(g)  Determine  an  equivalent  wave  height,  Hg,  for  the  total  fetch 
and  the  given  windspeed  and  water  depth. 


(h)  Calculate  the  fractional  growth  by 


Gi 


H 


am 


(18) 


(i)  Calculate  the  decayed  wave  height  at  the  end  of  the  fetch  by 


Hp  = Hm  - Gi  aVr  - H«^) 


(19) 
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2.  Wave  Period. 

As  waves  decay  over  the  fetch  segment,  the  significant  wave  period 
also  changes.  Very  long  waves  decay  rapidly;  shorter  waves  may  decay 
very  little  (see  Figs.  16  and  18).  This  means  that  the  significant  wave 
period  may  be  reduced.  As  a conservative  estimate,  it  will  be  assumed 
that  the  wave  period  remains  constant.  This  is  a conservative  estimate 
since  longer  period  waves  would  produce  higher  runtrp  on  a structure,  all 
other  variables  being  the  same. 

***************  example  problem  2******’******* 

GIVEN:  A coastal  area  is  flooded  by  a storm  surge  so  that  the  water 
depth  over  the  area  is  10  feet  (3.05  meters).  The  actual  fetch  across 
the  area,  in  the  direction  of  wave  travel,  is  3,000  feet  (914  meters). 
The  area  is  covered  with  thick  stands  of  tall  grass  and  a small  to 
moderate  amount  of  brush  or  low  bushy  trees  in  an  even  distribution. 


I 


I 
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The  windspeed  is  90  miles  per  hour  (132  feet  or  40.2  meters  per  second) 
and  the  initial  wave  height  at  the  seaward  edge  of  the  area  of  6 feet 
(1.83  meters);  the  wave  period  is  4.5  seconds, 

FIND:  The  decayed  wave  height  at  the  end  of  the  fetch. 

SOLUTION:  From  the  long  dash  line  in  Figure  1,  for  the  windspeed  of  90 
miles  per  hour  and  the  water  depth  of  10  feet 


ed  32.2  X 10 


= 0.0185 


giving  (at  the  intersection  of  the  above  line  with  the  long  dash  line) 
•SS  = 0.0075 

so  that  the  maximum  significant  wave  height 

= 0,.0075  (13211  = 4.1  f^et  . 

From  equation  (15), 

Hm  = 0.78d  = 0.78  (10)  = 7.8  feet 

and  from  equation  (16),  the  fractional  reduction  is 

D.  = ^ ~ = 7.8  - 6 _ = 0.486  . 

^ 7.8  - 4.1 

From  equation  (17),  the  equivalent  initial  wave  height 

ttsm  = 0-486  x 4.1  = 1.99  feet  ; 


from  Figure  1,  for 


M = 32 Cl-OPJ.  , 0.00368 
(132)2 


^ = 0.0185 
U2 


the  fetch  is  given  by 


; 


F = 760  feet  for  the  90-inile  per  hour  windspeed,  so  that  the  equiva- 
lent fetch  is 


Fg  = 760  feet 

The  vegetation  does  not  match  any  of  the  curves  in  Figure  15,  but  falls 
between  curves  B and  C.  Assuming  that  a moderate  amount  of  brush  will 
give  a friction  effect  about  halfway  between  the  two  curves,  from  curve 
B,  where  d = 10  feet,  fj>  = 0.20,  and  from  curve  C,  where  d = 10  feet, 
ff  = 0.485.  The  bottom  friction  is  then  taken,  in  this  case,  as  the 
average  of  the  two  values 


£ _ 0.20  + 0.485  _ n TA1 

2 ■ 


For  ff  - 0.01, 


ff  Ax 

0.01  X 6 X 3,000 

d2 

102 

for  ff  = 0.343, 

ff  Ax 

0.343  X 6 X 3,000 

d^ 

102 

for  T = 4.5  seconds  and  d 

= 10  feet. 

2ird 

_ 2j  (10)  _ 0.096 

gT2 

g (4.5)2 

From  Figure  16, 

Kj?. 01  * 0.80  for  ff  = 0.01  and  ff 

Kf^  = 0.105 

for  ff  = 0.343  and  f 

From  Equation  (13), 

1 - 0.105  0.895 

1 - K/.oi 

1 - 0.80  " 0.20 

from  equation  (12), 

II 

Q 

> 

X 

(1 

48  (3,000)  = 13,440 

(i.e.,  the  wave  decay  over 

3,000  feet  of  tall 

4.48  ; 


equal  to  the  wave  decay  over  13,440  feet  of  a sand  bottom  for  this  water 
depth  and  windspeed) . 


F = Fg  + F^  = 760  + 13,440  = 14,200  feet  . 
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For  a wind.speed  of  90  miles  per  hour  and  a fetch  of  14,200  feet  (from 
Fig.  1) 


^ = 0.0185  (as  previously  determined) 


F _ 32.2  X 14,200 


4 

l|2 


U 


(132)2 


= 26.24 


giving 


4 = 0.0071  . 

From  which  the  equivalent  wave  height. 


He  = 


_ 0.0071  U2  0.0071  (132)2 


32.2 


3.84  feet 


From  equation  (18),  the  fractional  growth  is 

^ He  3.84 

Gv  = = = 0.937  . 

'U  4.1 

The  decayed  wave  heiglit  is  then  given  by  equation  (19)  as 

(H^  - Hg^)  = 7.8  - 0.937  (7.8  - 4.  1)  = 4.33  feet  . 

At  the  end  of  the  fetch  segment,  the  wave  height  and  period  are  approx- 
imated by 

H^j  = 4. 33  feet 
T = 4.5  seconds 

it***1t****'k*******ir*it*it**ir*****it*‘it*it** 


VI.  SUMMARY  AND  CONCLUSIONS 

Tile  method  presented  in  this  report  gives  a first  approximation  for 
estimating  wave  heights  at  the  end  of  a fetch  with  a high  value  of  bottom 
friction  (e.g. , a flooded  area  with  dense  stands  of  grass  or  brush).  Only 
limited  data  are  available  for  wave  height  growth  or  reduction  for  waves 
passing  over  areas  with  dense  bottom  vegetation.  The  method  has  not  been 
verified. 


A substantial  amount  of  data  is  needed  for  waves  passing  over  areas 
of  flooded  vegetation.  The  method  should  be  verified  and  modified  as 
required.  The  shoaling  coefficients  (Fig.  16)  and  the  friction  factors 
(Fig.  15)  should  also  be  compared  with  values  from  actual  measurements. 
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